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Abstract: In this article a first step is made towards the extension of Britto-Cachazo- 
Feng-Witten (BCFW) tree level on-shell recursion relations to integrands and integrals 
of scattering amplitudes to arbitrary loop order. Surprisingly, it is shown that the large 
BCFW shift limit of the integrands has the same structure as the corresponding tree 
level amplitude in any minimally coupled Yang-Mills theory in four or more dimensions. 
This implies that these integrands can be reconstructed from a subset of their 'single 
cuts'. The main tool is powercounting Feynman graphs in a special lightcone gauge choice 
employed earlier at tree level by Arkani-Hamed and Kaplan. The relation between shifts of 
integrands and shifts of its integrals is investigated explicitly at one loop. Two particular 
sources of discrepancy between the integral and integrand are identified related to UV and 
IR divergences. This is cross-checked with known results for helicity equal amplitudes at 
one loop. The nature of the on-shell residue at each of the single-cut singularities of the 
integrand is commented upon. Several natural conjectures and opportunities for further 
research present themselves. 



Keywords: Amplitudes. 



Contents 



1. Introduction 1 

2. BCFW recursion for tree level amplitudes 4 

2.1 Powercounting Feynman graphs in a special gauge 6 

2.2 Adding scalar and quark matter 10 

2.3 Remarks 11 

3. BCFW shifts of loop amplitude integrands 12 

3.1 String theory perspective 13 

3.2 Shifts of the integrand at arbitrary loop level 14 

3.3 Remarks 16 

4. Difference between integrand and integral at one loop 17 

4.1 BCFW shift of the integrals at one loop 17 

4.2 One loop momentum integrals contributing to large BCFW shift 21 

4.3 One scalar loop 24 

4.4 Remarks 27 

5. Towards BCFW on-shell recursion relations at loop level 28 

6. Discussion and conclusions 31 

A. BCFW shifts of fermions in D dimensions 32 

B. One loop diagrams contributing to the BCFW shift at leading order 37 



1. Introduction 

Scattering amplitudes are on the crossroads of many developments in physics and are 
therefore prime objectives for the developments of new calculational methods. In recent 
years the technology for analytic calculation of these amplitudes has seen many advances 
inspired by Witten's twistor insights [1]. Contact of these developments with the calculation 
of experimentally (still!) relevant NLO calculations is being made. However, many of the 
recently developed techniques are applicable for QCD only at the tree or one-loop level, 
while many of the most exciting analytic developments have been at high loop level in 
J\f = 4 super Yang-Mills theory. NNLO calculations in QCD are certainly on the (longer 
term) wishlist of experimenters. Apart from this direct motivation from the experimental 
side there is the inherent intuition that many of the recently found hidden structures 
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in maximally supersymmetric Yang-Mills might have some form of non-supersymmetric 
counterparts. In this article the full extension of one particular tree level technique to the 
loop level at in principle arbitrary loop order will be initiated. This technique will be that of 
the Britto-Cachazo-Feng-Witten (BCFW) on-shell recursion relations [2], [3]. In a nutshell, 
at tree level these relations allow one to calculate any amplitude from the analytic form of 
the three point amplitude only. There is a broad parallel here to bootstrap equations in 
CFT as pointed out in [4]. This looks especially natural from a twistor transformed point 
of view [5] , [6] . From the CFT bootstrap perspective it is natural to widen the investigation 
of recursion relations to the loop level. 

On-shell recursion at loop level for complete amplitudes has certainly been discussed 
previously in several places in the literature, starting in [7], where the purely rational 1- 
loop amplitudes in pure Yang-Mills were studied starting from their known expressions. 
In further work this has been developed into a method for calculating rational terms [8], 
[9] in general which has in addition been implemented in numerical code [10]. Recursion 
for coefficients of the loop integrals has also been studied [11] and a recent application 
of these techniques includes for instance the calculation of rational terms in some pure 
Einstein one loop gravity amplitudes [12]. One technical problem in all these approaches 
is in general the appearance of so-called boundary contributions for large BCFW shifts 
which spoil direct predictability. Controlling these can be done using auxiliary recursion 
relations [8], [9] and general powerful consistency conditions but a direct understanding of 
the boundary contributions has been missing so far. Boundary contributions appear also in 
other context, see e.g. [13], [14] for a very recent discussion. Moreover all of the mentioned 
work depends on properties of amplitudes special to one loop. 

In this article the known issues with on-shell recursion for amplitudes at loop level are 
first sidestepped by considering the integrand of the amplitudes instead of the integrated 
expressions. It is quite straightforward to check in several explicitly known examples of 
amplitudes in gauge theories that the scaling of the integrand matches that of the corre- 
sponding tree level amplitudes. This can be seen in for instance the case of MHV ampli- 
tudes in N = 4 at one loop as obtained by long ago [15] and the all-plus amplitudes in 
pure Yang-Mills, also at one loop (using the fact that these amplitudes are related [16]). 
By inspection the full integrand of these amplitudes is seen to behave under BCFW shifts 
exactly as the tree amplitude. The same result also follows for the integrals appearing in 
the planar 4 point amplitude in N = 4 at five loops as given in [17]. This leads to the 
following suspicion: 

Suspicion 1.1 In four dimensions and up, as long as all loop momenta and the external 
momenta are sufficiently generic, the integrand of Yang-Mills and Einstein gravity ampli- 
tudes coupled to fundamental or adjoint matter scales under a BCFW shift of a pair of 
outside legs as the corresponding tree amplitudes to all loop orders. 

Although the focus in this article will be on gauge theory, Einstein gravity is included 
in the above suspicion since this is quite natural to suspect. For gauge theory there is one 
class of shifts which does not exist at tree level: those of gluons on different color traces. 
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One expects these to behave as non-adjacent shifts at tree level, i.e. they should be more 
(at least one power of l/z) suppressed under large shifts. As explained below the shift 
behavior proves the integrand can be reconstructed from the propagator pole singularities 
in the shift parameter. If the residues at the poles of the integrand can be interpreted in 
terms of amplitudes again this implies 

Suspicion 1.2 In four dimensions and up, as long as all loop momenta and the external 
momenta are sufficiently generic, the integrand of Yang-Mills and Einstein gravity am- 
plitudes coupled to fundamental or adjoint matter obeys (a form of) on-shell recursion 
relations to all loop orders. 

The residues at the poles in the integrand correspond to single cuts: replacing one 
particular propagator by a delta function. These single cuts of the integrand are reminiscent 
of but distinct from Feynman's tree theorem [18] or even more the streamlined versions of 
the tree theorem in [19], [20] and especially [21]. The difference is that the construction 
studied in this paper will turn out to involve a certain strict subset of single cuts only 
while keeping propagators the same as usual. In contrast, the mentioned applications of 
Feynman's tree theorem study all cuts in the original form or all single cuts in the more 
advanced version with certain modified propagators (or boundary conditions). Note that 
for direct numerical applications the approach of Catani et. al. might be better suited as it 
treats all singularities equally, allowing in principle to treat all the IR singularities of real 
and virtual corrections on an even footing. The usefulness of applying a more conventional 
expand-in-a-basis unitarity-type approach using single particle cuts has been studied in 
[22]. 

From the outset it is clear that a tension exists between the limit of large BCFW 
shift of the integrands versus that of the integrals. This arises since the integrals are only 
well-defined after (dimensional) regularization, so this can be viewed as an order of limits 
problem. For the leading poles in the dimensional regularization parameter e this is not 
expected to be much of a problem, but for sub- leading terms a problem might arise. At 
one loop this is a well-known phenomenon related to the appearance of so-called 'rational 
terms'. This is motivation to study the difference between integrands and integrals more 
closely as will be done below. 

This article is structured as follows: in section 2 the standard tree level BCFW recur- 
sion relations are first reviewed. For the derivation the large BCFW shift behavior of the 
tree level amplitude is crucial and a completely diagrammatic version of this derivation is 
presented using powercounting in a special lightcone gauge. This argument is independent 
of much of the fine details of the exact correlation function under study such as the loop 
order. This is expanded on further in section 3. In particular, the diagrammatic argument 
shows that the integrand of Yang-Mills theory coupled to various kinds of matter has the 
same BCFW shift as the tree level amplitude for shifts of particles on the same color trace 
thereby establishing suspicion 1.1 for gauge theory for this class of shifts. This result on 
integrands begs the question how the shift relates to the integrals of the integrands. This 
question is studied first at one loop in section 4 using the background field method to sep- 
arate gauge choices of trees and loops. In particular the large BCFW shifts of the purely 
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rational amplitudes at one loop are studied. In section 5 a preliminary discussion of the 
relation of the single cuts of the integrand to lower point or lower loop amplitudes is given 
as a first step in the direction towards suspicion 1.2. A general discussion and conclusion 
section rounds off the presentation. In appendix A shifts of other particles than gluons are 
considered in four and higher dimensions. Appendix B contains a set of explicit results for 
a series of Feynman graphs used in section 5. 

Note added in proof: Simultaneously with this article closely related and independent 
work by Arkani-Hamed et.al. [23] has appeared on the archive which contains an explicit 
example of a recursion relation for the integrand for planar N = 4 super Yang-Mills theory 
in four dimensions. 



2. BCFW recursion for tree level amplitudes 

As BCFW observed amplitudes can be turned into a function of a single complex variable 
by shifting the momenta of two legs of the amplitude as 

hi V k% + qz /„ 

rvj r K/j QZ , 

while requiring the vector q to obey 

q-ki = q-kj=q-q = . (2-2) 

which keeps the invariant lengths kf and fc| invariant. These equations have two complex 
valued solutions for q, as can be checked in an appropriate ('center of mass') Lorentz- 
frame. A covariant solution can also be constructed using higher dimensional spinor helicity 
methods [24] [25], see appendix A for details. The shift turns the amplitude into a function 
of the complex shift parameter z. This function can then be contour integrated to yield 
the amplitude of interest, 

A(0) =1 . (2.3) 

Jz=0 z 

Up to this point no assumption about the loop order of the amplitude has been made 
up to the reasonable assumption of the singularity at z = to be a single simple pole. If 
the amplitude is tree level however, it must be a rational function of the momenta and 
hence also of z. In this case one can pull the contour to infinity as depicted in figure 1 to 
yield 1 

A(0) = <f ^ =J2(z = finite) + ^ (z = infinite) . (2.4) 

J 2=0 Z 



res 



In this formula the finite z poles have known residues by tree level unitarity. These residues 
are a product of tree amplitudes, summed over all states in the cut channel and summed 
over all factorization channels where the shifted particles are on distinct tree amplitudes 
in the residue (possibly taking into account color ordering). Importantly, the involved 
amplitudes have a strictly smaller number of particles. The infinite z residues do not have 
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Amplitude 




Figure 1: Sketch of contours and singularities used for BCFW recursion at tree level in equation 
2.4 
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Table 1: leading asymptotic power in z K of the adjacent BCFW shift of two gluons in a tree 
amplitude for all possible polarizations of the shifted gluons 

a similar interpretation. If the latter residues are absent however, equation (2.4) is an 
on-shell recursion relation. 

The challenge is therefore to derive the large z-behavior to show the absence of these 
residues. In general this depends on the helicity of the particles being shifted and whether 
or not the shifted particles are adjacent in the color ordering (for gauge theories). For 
gluons in Yang-Mills theory the large z-behavior is for instance displayed in table 1 for 
shifts of particles adjacent in the color ordering. In this table the difference between T and 
T2 is whether the inner product between the two polarization vector vanishes or not. If it 
does, the shift is better behaved. Non-color-adjacent shifts are suppressed by an additional 
factor of (-). See table 2 in particular and appendix A in general for how amplitudes 
behave for shifts of gluonic and fermionic legs. . 

Table 1 can be derived through various means but perhaps most physical Arkani- 
Hamed and Kaplan [26] pointed out that a large BCFW shift basically describes a hard 
particle shooting through a soft background given by the other gluons (see equation 2.1). 
Technically, this can be described handily by the background field method. The background 
field can be put in the q lightcone gauge to show that the naive powercounting of the large 
z-behavior of diagrams in which the shifted legs are on different vertices is suppressed by 
This leaves local vertices which can be analyzed by expanding the Yang-Mills action 
coupled to for instance fermions to second order in the 'hard' fields using this method gives 

£ = D u a^D u a^ + a^a u F^[A] + a M (^ 7 ^ + VV^) + A^fij) , (2.5) 
1 Here and in the following the § symbol stands for the integral ^— § dz. 
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where and tp describe the hard particles, while A^ and \& describe the background field. 
D u is the covariant background field derivative. The linear expansion term, 

" C[A^])+^(4rC[A^}) , (2.6) 



vanishes by the equation of motion for the soft fields. The structure of this Lagrangian 
immediately leads to a prediction of the large z structure of the two particle current which 
in turn yields table 1 when paired with the external wave functions. 

2.1 Powercounting Feynman graphs in a special gauge 

The original AHK argument depends on an application of the tree level equation of motion 
for the external fields. Below we show one can reproduce the same argument in a more 
diagrammatic language, which will be useful below. The BCFW on-shell recursion relations 
have been studied from the diagrammatic point of view before in [27], but there the focus 
was on the much harder task of understanding the terms in the recursion relations directly. 
Here the focus will be firmly on the large BCFW limit. As a bonus explicit expressions 
for the first and second non-vanishing coefficients in the large z limit will be produced. 
For this the amplitude under study is computed in what will be referred to as the Arkani- 
Hamed-Kaplan (AHK) gauge, 

= . (2.7) 
The corresponding lightcone gauge propagator reads 



G{k) 



1 



k ■ q 



(21 



with 7] the flat space metric. This propagator connects the usual (color ordered if required) 
three and four point vertices and the external fields. The auxilliary field is left in the 
Lagrangian in contradistinction to the usual lightcone approach. For our purposes the 
external wave functions of the shifted legs will be kept in a different gauge than the other 
legs which will be put in q lightcone gauge. This is necessary since in the large z limit 
the momenta of the shifted legs become proportional to q and the AHK gauge is therefore 
singular for these two special external wave functions. 

There is a more serious gauge singularity for those diagrams for which the shifted legs 
end on the same three vertex. For this special class q is orthogonal to the momentum in 
the soft leg, 

q ■ (fci + fej) = , (2.9) 

and the above lightcone gauge propagator in equation (2.8) is singular. This difficulty can 
be circumvented by first imposing the auxilliary lightcone gauge 

(9 + 3*1)^ = 0. (2.10) 

For the diagram which is singular in the limit x — > the relevant propagator shows of 
course an explicit pole in x, 

G{kl + k2) ^ = (hTW{^ xjk^) J ' (2 - n) 
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This pole should be a gauge artefact since the complete amplitude is gauge invariant. This 
can be made more precise. Prom the usual expression of the (color ordered) three vertex in 
Yang-Mills the contribution of the diagrams with the shifted legs ending on the same three 
vertex coupled to an arbitrary remainder through the above lightcone propagator can be 
written as 

~ e%<% [ V ^(2zq + h - k 2 ) p + + k 2 ) u - jf v {k x + G K p 



6i e 



l c 2 



(fci + k 2 



\2 



7]^(2zq - 2k x ) K + (rfi(ki + k 2 ) u - rfcfa + k 2 )^)+ 



[{q)^(h + k 2 )y ~(q + xk x )u{kx + k 2 ) fi ] (fci + k 2 ) K ) 



(2.12) 



x(k\ ■ k 2 ) 

where the hats on the momenta indicate BCFW shifted quantities. To obtain this one uses 

7k K 

{kl -k 2 YG{k 1 + k 2 T p = - w ^ w . (2.14) 

All other (parts of the) diagrams in the theory do not have poles in x for sufficiently 
generic momenta. Hence one expects the poles on the last line of equation (2.12) to vanish. 
Actually, since the possible pole arises from either 

<7M e i or <l» e 2 . ( 2 - 15 ) 
it is seen that the possible pole in x is lower order in z since 

q fl ef = ± i J ^ef i = l,2, (2.16) 

where the sign depends on the particle. Hence for our purposes the terms involving ^ 
poles can be consistently dropped. Moreover, the entire pole at x = vanishes by gauge 
invariance of the one particle off-shell current which contracts into the above vertex. This 
follows as the momentum of the off-shell leg is {k\ + ^2)- 

Hence the class of diagrams where the shifted legs end on the same vertex contributes 
terms of the form 

A(z) -> e v 2 (k 2 ) ilV = 4 (z) e v 2 (z) (z^h UA) + 0-M +o(^\\+... , 

(2.17) 

to the large shift limit. Here ji^v is an antisymmetric matrix and the dots indicate 
the diagrams where the shifted legs end on different vertices to be analyzed below. The 
structure of the leading and sub-leading terms here is basically the AHK result [26]. The 
functions /j are in general polynomials in with non-trivial constant term. 

For the class of diagrams just considered the function f\ for instance can be calculated 
from the coupling of the one leg off-shell current to the above three vertex current with 
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two shifted legs as 

fl = U - ^) JftX) ■ (2-18) 

Here is the current and X stands for the quantum numbers on the other legs of the 
diagram. Note that in a real sense the leading-in-z contribution to the large BCFW shift 
can be thought of as an effective particle made out of particles one and two. This effective 
particle interpretation follows also from the CFT analysis of BCFW shifts in string theory 
[28] [4]. In the string context the corresponding particle is sometimes referred to as the 
pomeron [29]. 

For the above result the contributions to the amplitude of those diagrams where the 
shifted legs end on different vertices has been disregarded. This will now be justified. The 
z dependence of these diagrams can be obtained from a refined form of powercounting. 
Note that in every tree diagram there is a unique path from one shifted leg to the other. 
This path will be referred to as the 'hard line' in the following as it is the path along which 
the z dependence flows. Every three vertex along the hard line scales as (z 1 ) with the 
legs contracted into q and a metric while every four vertex scales as with two metrics 
contracting the legs. Along the hard line the lightcone gauge propagator scales as 



G(k) 



1 

¥ 



k ■ q 

1 



k 2 ± 2zqk 



Qfikv + k^q v 
k ■ q 



2z 

=F 



k 2 ± 2zqk 



k ■ q 



(2.19) 



which contains a term proportional to z° in the limit. The sign depends on the routing of 
the momentum. However, although the superficial degree of scaling of a given diagram is 
now z" 3 with v% the number of three vertices, this is never realized. The above propagator 
remains after all orthogonal to q, so the part of the three vertex proportional to q never 
contributes. The order z part of the above propagator has to be contracted with soft field 
propagators which are orthogonal to q or into the soft external fields also orthogonal to q. 
This leaves the two hard external fields and contractions along the hard line. The only way 
a possible z-dependence could arise from the external hard wave functions would entail 
terms like 

e1(h)kl = eKk) {% + k]) = ef(fci) (*J + $ ■ (2-20) 

and its natural conjugate which are therefore suppressed. Moreover, for this class of con- 
tributions there is now an external field contracted into the lightcone gauge propagator 
along the hard line which leads to a (j\ suppression. It is seen that only contracting the 
momenta on two three vertices using the lightcone propagator and with the shifted fields 
on the external legs of these three vertices yields an order z° contribution. In fact, the 
form of this is easy to derive up to a constant, 

4(k)e v 2 {k)V^ K G KX Vl a ~ z°e1e^ KVva + O 0) . (2.21) 

Hence this can be treated as an effective 4 vertex. At this order in z there is also a direct 
contribution of the usual Yang-Mills four vertex. These two terms can be summed to read 
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in the large z limit in color ordered form 

4 [k) e v 2 (k 2 ) V^ pa ~ (kx) e v 2 (k) (v w Vua ~ M^) , ( 2 - 22 ) 

where a and p connect to the rest of the diagram. Note that this is anti-symmetric in the 
indices /i and v (as well as /x and a), although the two four vertex terms separately are not 
anti-symmetric 2 . 

Summarizing, in the limit of large BCFW shift of two color adjacent gluons the scat- 
tering amplitude scales as 

A(z) -> e£ (%) e\ (fc 2 ) AV = e^ (z) e v 2 (z) (z^ iv h{l/z) + f 2>flu (l/z) + O Q 
with /i given in equation (2.18) and the antisymmetric term f 2 given by 



(2.23) 



2,/j.u 



(VfJUtikl + k 2 ) v - T) VK (ki + k 2 )^)Jg + ^ i^tMpVva ~ VupV/ia) J q{ X i) J q{ X i+l, 



(2.24) 

where J q are again the currents calculated in AHK gauge and the sum runs over the ways 
to subdivide the set X of external particles into two subsets Xi,Xi + i, keeping the color 
order of the particles. The structure displayed in equation (2.23) is equivalent to the result 
of [26] but bypasses the use of the tree level field equations. Using the explicit scaling 
of the external wave functions reviewed in the appendix the analysis just presented leads 
directly to table 1. 

Comparison to Feynman-'t Hooft gauge 

The above large z contributions although obtained in a particular gauge, are expected to 
be gauge invariant since are the limit of a gauge invariant quantity. As a cross-check the 
leading order in z part of the above result also follows from power-counting the diagrams 
which contribute at this order in Feynman-'t Hooft gauge. What changes here is the form 
of the propagator as well as the conditions on the fields connecting to the hard line. The 
leading z contribution for Yang-Mills theory in this gauge is again formed by a hard line 
which only consists of three vertices. By the general structure of the three vertex the hard 
line will have two fields external to the hard line contracted with a metric, while the other 
fields are contracted into q p : allowing other fields will lower the power of z. Therefore, for 
pure gauge theory the hard line contributes 

~ [z]A 1 ^ + Kyfh, + qjt v + fiqtfj) A 2 >» , (2.25) 

in the large z limit, where and A 2 connect to the external wave functions. The quantities 
K and K are arbitrary vectors which stand for further connections in the diagram. From 
equation (2.16) it follows that the leading contribution from any diagram with a hard 



2 It is often said that BCFW recursion relations make the four point vertex of the Yang-Mills Lagrangian 
obsolete. While strictly true, the computation here shows clearly that in deriving the relations the four 
vertex is crucial. 
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line is proportional to the metric contraction in pure Yang-Mills, while the proportionality 
factor is calculated in the same way as derived above in AHK gauge by contracting with 
the current. The powercounting argument in Feynman-'t Hooft gauge could probably be 
extended to the subleading contributions at order z° with more effort. 

2.2 Adding scalar and quark matter 

For theories with more general matter content than tree level Yang-Mills matter couplings 
are important. The effects of these on shifts of gluons can be analyzed by a simple extension 
of the above. 

Scalar s 

Scalar matter will be important since ghosts for instance generically have scalar type cou- 
plings to glue. With the gluons in AHK gauge there is a limited number of new diagrams 
appearing for the BCFW shift of two gluons if coupling to scalar matter is included. These 
involve one or more scalar hard lines. For every connected scalar hard line it is clear that 
this part of the hard line will contribute to leading order 

~A^A v (zq»q v + K»q v + K v q» + o(^X\ , (2.26) 

where the explicit gluon fields are simply placeholders for any connection further along the 
hard line and K and K are again arbitrary vectors which stand for further connections in 
the diagram. In AHK gauge it should be noted that these vectors would always involve 
either a q^ connecting to a soft or to a hard line. For soft lines these diagrams actually 
vanish while for hard lines this contribution is suppressed by (-). From the analysis of the 
gluon diagrams in AHK gauge it is clear that only those diagrams where the scalar line 
contains just zero or one element and where the outside legs connect directly to the scalar 
line will potentially contribute to the large shift limit. The diagrams with one element 
would always contain at least one of 

k],4 or ky 2 , (2.27) 

as terms which have a possible z dependence in the denominator, but these terms vanish 
of course. The only remaining contribution is formed by the two hard lines connecting 
directly to the 2 gluon 2 scalar vertex. This is of course proportional to the metric and of 
order (z )- The effects of this vertex should be included in equation (2.23) if one wants to 
calculate the exact coefficients of the large shift, with of course appropriate scalar currents 
attached. 

With the gluons in Feynman-'t Hooft gauge instead of AHK gauge equation (2.26) 
continues to hold. This does not change the form of the current in equation (2.25) when 
scalars are included and hence therefore does not change the form of (2.23). 

Quarks 

Fermionic spin one half contributions to the hard gluon line are suppressed by one power of 
z compared to the leading contribution in the gluonic case. This follows because the powers 
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of z in the numerator are generated by the fermionic propagators, not by the vertices. In 
the case of minimally coupled spin one half matter, there are only three particle couplings. 
These lead for the diagrams which form a fermionic path along the hard line to 

~ ApAjty f 7m G7„ + O f^j V , (2.28) 

where ip and £ are arbitrary spinors which symbolize the coupling to the soft parts of the 
diagrams and the matrix G reads 

I 

G = U)U (7« «)) > (2-29) 

j=Q 

where the indices fij connect to soft parts of the diagram and there are I external soft fields 
attached to the fermion line. Since (j^) 2 = one can also write this as 

G=d)\{{q H ) ■ (2.30) 

3 

The resulting gamma matrix structure can easily be split into a symmetric and anti- 
symmetric part to give 

where square brackets denote antisymmetrisation. In AHK gauge any non-trivial con- 
tract to the soft fields and the diagrams vanish, leaving the contribution with a single 
element. Just as above it is easy to see the hard fields have to attach directly to the 
fermion line with just one propagator to have a non-trivial contribution. This allows to 
drop the second symmetric term in the above as this is subleading in z after contraction 
with the external hard wave functions. The remaining two terms must be added to (2.23) 
to calculate exact limits. 

With the gluons in Feynman-'t Hooft gauge it is easy to see that just as for scalars 
spin one half matter doesn't change the overall structure of the leading large z limit on 
two gluons. 

2.3 Remarks 

Adding scalar fermion couplings does not change the preceding analysis as this coupling is 
automatically suppressed. It should be stressed that the derivation of the leading z behavior 
as presented in this section applies to basically any minimally coupled gauge theory with 
scalar or fermionic matter. Also the above clearly shows that at leading order it is a 
property of a subset of the full Feynman graphs: only the hard line needs to be considered 
with a simplified set of Feynman rules. This shows that on-shell recursion techniques can 
be applied to more general correlation functions than those used for scattering amplitudes 
only. Also, since the leading large z behavior arises from a either a three vertex diagram 
with the shifted legs attached directly or from a diagram with the shifted legs connected 
by one propagator only, it is easy to see that non- adjacent shifts are suppressed by at least 
an additional (-). 
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Bonus relations? 

Note that knowledge of the exact form of the large BCFW shift limit of the scattering 
amplitudes can in principle be used to derive new relations between scattering amplitudes 
at tree level. This follows because the BCFW derivation of equation (2.4) can now be 
modified as 

0=6 (l + az)A(z) = ^(z = finite) + ^ (z = infinite) , (2.32) 

z — res res 

for some arbitrary coefficient a. Here the modified finite z residues are now a sum over 
products of tree amplitudes times a momentum dependent function. This sum equals the 
now non-vanishing residue at infinity which can be calculated explicitly for the (+— ) good 
shift through the above diagrammatic techniques for instance. For shifts of a pair of helicity 
equal gluons there is a gap consisting of the missing explicit order Q) contributions from 
the diagrammatic analysis for this particular shift. It would be interesting to explore this 
further as the currents can be calculated using auxiliary recursion relations if needed. 

Currents 

In the above one leg off-shell currents appear, calculated in two different gauges. However, 
it is fairly easy to see that the currents themselves are largely gauge invariant as they can 
be calculated through an auxiliary BCFW recursion relation. For this one simply uses an 
auxiliary BCFW shift. The resulting recursion relations, when iterated, only rely on a new 
three vertex: the three particle current with one leg off-shell. Since for this vertex 

(KYJl = (2.33) 

holds with K the momentum in the off-shell leg, it holds for the n-particle current as well. 
Using off-shell recursion relations this was proven in [30]. Note also that the current is 
invariant under gauge transformations on each leg separately, i.e. the current is invariant 
under 

6 -»• & + h . (2.34) 

The same conclusion follows by a short modification of the argument above for the singular 
diagrams in the AHK gauge: requiring absence of gauge irregularities in the x — > limit 
shows that the current must be gauge invariant if the parent amplitude is. 

3. BCFW shifts of loop amplitude integrands 

In this section the discussion is extended to loop level. There are some obstructions to 
applying the above tree level derivation of on-shell recursion relations and in particular of 
equation (2.4) to the loop level directly. Most obviously, at loop level there will be branch 
cuts in the z plane in general such as those sketched in figure 2. 

Note however that the contour integral at z = can still be deformed to encircle 
the branch cuts. Assuming the branch cuts do not overlap, one can see that the integral 
around the branch cuts is simply an integral over the imaginary part of the discontinuity 
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across the branch cut. Since the discontinuities in the branch cuts in amplitudes generically 
correspond to cut propagators in the Feynman graphs, one sees immediately that there is 
scope for a loop level analogue of on-shell recursion. For this to work properly, the integrals 
over the arc at infinity must still vanish or at least give a well-defined result. Hence for 
immediate applications one still would like the amplitude to scale well at z = oo, as long 
as z is not directed along one of the branch cuts. 

Sidestepping integrated amplitudes first, in this section the previous discussion is ex- 
tended to loop level integrands. This is natural as the integrand is just as the tree level 
amplitude a rational function of the momenta. Hence it scales as a rational function of 
z. The main argument is based on powercounting, but the main idea of the present sec- 
tion first arose from a simple consideration in string theory which will be presented first. 
Readers not interested in string theory should probably skip the next subsection. 

3.1 String theory perspective 

On-shell recursion at string tree level has recently been proven using a certain limit of the 
OPE between two adjacent vertex operators [28] [4]. This argument can be applied to a 
general string theory amplitude which reads for the closed string for instance 

A n = j dm{V x {zT L )...V n {z n )) , (3.1) 

where d[i is the measure on the moduli space of spheres with I handles and n operator 
insertions. Despite its complicated form the moduli space integral does not depend on the 
momenta of the particles. 

On a genus g Riemann surface the OPE between for instance two closed string tachyon 
vertex operators in the bosonic string reads 

:V(zi)::V(z 2 ):~ e G «( fc i'*»>(*-*0 :V(z 1 )V(z 2 ): , (3.2) 

where G g is the solution to the equation 

d w dyjG g (w) = 5 2 (w) (3.3) 

on the Riemann surface of interest. Crucially, for small w G g tends to the tree level two 
point correlator on the sphere, 

G g {w) ~ log(io) + log(W) + 0(w) . (3.4) 
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It is easy to see that the prefactor on the right hand side in equation (3.2) is shift inde- 
pendent because it involves k\ ■ &2- I n other words, if we assume just as at tree level that 
wz ~ 1 and take the large z limit, the analysis quickly reduces to the analysis of BCFW 
at string tree level discussed in [28] and [4]. For more general shifted particles in both 
open and closed (super)strings a similar derivation also holds. Ignoring the measure and 
possible complications from singular points on the moduli space, the above suggests the 
natural extension of suspicion 1.1 to string theory: 

Suspicion 3.1 In four dimensions and up, away from singular points on the moduli space, 
the integrand of open or closed string theory amplitudes in a flat background scales under 
a BCFW shift of a pair of outside legs (at least as good) as the tree string amplitude to all 
loop orders. 

Taking the naive field theory limit of this string theory suspicion gives the above 
field theory suspicion (1.1). This explicitly ignores any possible complications with the 
decoupling of massive modes for instance. Note that the above argument holds for adjacent 
vertex operators only in the case of the open string. Non-adjacent shifts would be expected 
to be exponentially suppressed. 

3.2 Shifts of the integrand at arbitrary loop level 

Note that even before beginning it is expected that the integrand of loop level amplitudes 
will scale as a rational function of z. This simply follows from any representation of the 
integrand of the loop amplitude in terms of Feynman diagrams. The large z shifts will 
involve shifts of the loop propagators by 



{li + h + K) 2 + m? (h + h + K) 2 + 2zq ■ (h + K) + m? ' y ' 1 

for a shift on momentum k\ within a loop made of a particle of mass nii and loop momentum 
li, where K is additional unshifted momentum (which explicitly rules out fo). As long as 
the loop momentum and the external momenta are real and generic 

q-(h+K)^0 (3.6) 

will hold and naive powercounting can be applied. Whenever the counting of z factors 
employed below will fail is where 

q-(li + K) = 0, (3.7) 

or where the loop momentum is infinite. For real momenta the condition (3.7) defines a 
D — 2 dimensional sub-manifold for the loop integration. As long as the integrand does not 
have delta-function support on this manifold (which it will not have for generic and real 
external momenta), these regions will not contribute to the amplitude. More generally, one 
can study points where the loop momentum is also sufficiently generic. This allows one to 
bypass the reality constriction on the loop momenta. 
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For the argument below to be transparent it is convenient to choose a convention 
for how the momenta in the loop depend on the external momenta: after all this can be 
changed by a linear, finite shift of any of the loop momenta. This shifts the integrand by 
a function which integrates to zero. For convenience one can choose the convention for the 
routing of the loop momenta such that the z-dependence in the loop diagram flows along 
a minimal path between the legs with shifted momenta, i.e. along the hard line. That 
such a routing exists for shifts of two particles within the same color trace can be seen by 
following the path on the color trace: the routing condition determines a convention for 
each loop encountered along the path. This fixes the z dependence of all loop momenta, 
leaving non-z-dependent loop momenta undetermined. 

For shifts of particles on different color traces no canonical choice of hard line exists. 
At one loop shifts on different color traces can be related to shifts of particles on the same 
color trace through the relations derived in [15] . Hence this is only an issue at two loops and 
beyond this case will not be studied in this article, save for the formulation of the suspicion 
that shifts of particles on different color traces are suppressed by at least one power of 
1/z compared to a shift of color adjacent particles. This depends on the ability to choose 
the loop momentum dependence on z to be only along a particular line in the diagram: 
following through the diagrammatic analysis this yields an 1/z suppression. Again, this 
case will not be further considered here. 

Powercounting 

The crux of the argument in this section is that the above diagrammatic analysis of the 
shift behavior in AHK gauge although motivated at tree level did not depend on on-shell 
constraints on the soft external legs at any point: it applies to a wide class of correlation 
functions. In particular, it does not depend on the loop order in an amplitude calculation. 
Up to a minor subtlety discussed below this will lead to the conclusion that the same 
analysis presented holds for the sum over diagrams in this particular gauge at any loop 
order, and in particular the integrand shifts according to the general structure outlined in 
equation (2.23). Here the functions f\ and ji in this Ansatz will not be given by (2.18) 
and (2.24) any more but by suitable generalizations. Hence the integrand of gauge theories 
coupled to various forms of matter scales as in table 1 to all loop orders for shifts of 
color adjacent gluons. The analysis of shifts of fermions considered in the appendix also 
generalizes immediately. 

There is a subtlety in the above reasoning which has to do with singular diagrams in 
the AHK gauge: there is an additional class of diagrams which consist of (generalized) 
'triangle' shaped loop diagrams on which the shifted legs end directly and only one leg 
connects to the rest of the diagrams. Just as at tree level this is singular in AHK gauge 
because the momentum k\ + &2 flows through the off-shell leg. In fact, this class has two 
different elements. One is formed by coupling the hard fields to the effective four vertex 
and tying the remaining two legs at the required loop order into a single external field. 
This diagram is of order (z°) and either proportional to the metric or anti-symmetric by 
the structure of the effective four vertex. 
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The other is formed by the tree level three vertex coupling to the hard fields and with 
a propagator correction graph on the external line. This has the same form as the tree 
level contribution (order (z 1 ) and (z )). This can be seen in a bit more detail. Since the 
propagator correction has to be transversal, the denominator of this correction must have 
the tensor structure 

Tftf - + kzYih + k 2 y , (3.8) 

to any loop order. Contracting with two propagators in the modified AHK gauge as above 
gives simply a tensor structure of the lightcone gauge propagator plus an extra piece of the 
form 

This extra term does not contribute since it contracts into the soft part of the diagrams. 
This reduces the analysis to terms studied before at tree level. 

Note that calculating loop integrals in a lightcone gauge directly can be fraught with 
difficulties as poles of the form 

~ • (3J0) 

arise which can however be overcome, see [31] and [32]. These poles cancel within the full 
amplitude expression as they should by gauge invariance, but can be hard to treat directly. 

In Feynman-'t Hooft gauge the leading order (z 1 ) behavior of the integrand can also 
be reproduced following the same reasoning as at tree level. Similarly, in this gauge the 
sub-leading behavior is harder to see. 

3.3 Remarks 

There are a number of natural expectations which follow from the above result. In brief, 
for every theory for which the BCFW shift behavior has been analyzed at tree level one 
expects the same behavior for the loop level integrand, to all loop orders. In eft 4 theory 
for instance this is easy to verify as the on-shell conditions for the external particles have 
no impact on the powercounting. For Einstein gravity the same is expected to hold, but 
this is although much more interesting also much harder to prove. As touched upon above 
as motivation, it is expected that for string theory a similar phenomenon holds for the 
integrand. 

Most of the explicit results above are for color- adjacent gluons. Note though that the 
argument does not depend on the color details of the rest of the diagrams: it holds for shifts 
of any two color adjacent gluons. In particular it holds for non-planar diagrams. Further, 
one can generalize easily the above to shifts of gluons on the same color trace as there still 
is a canonical hard line along the edge of the diagram. As discussed before, shifts of gluons 
on different color traces beyond two loops will be left to future work. 

A natural extension of the results of this section is to use coherent state methods 
(BCFW 'supershifts') in maximally supersymmetric theories in four dimensions [33], [34] 
or higher [25]. This can also be used for instance to derive a slate of predictions for shifts 
of any pair of members of the A = 4 multiplet, at any point of the moduli space [35]. 
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4. Difference between integrand and integral at one loop 

After the discussion of the BCFW scaling of the integrand it is a natural question how 
this extends to the scaling of the integrals of these integrands. After all, one would like to 
study recursion for the amplitudes directly as these are the physical objects. The difference 
between the two scalings depends on the ability to interchange integration and limit. This 
can be phrased as an order of limits problem: the loop integrals are regularized to make 
them well-defined by, say, dimensional regularization. Then the question is whether the 
limit of large BCFW shift (z — > oo) commutes with the limit in which the dimensional 
regularization parameter vanishes (e — >■ 0). For the leading singularities of the integrands 
this is (almost) bound to work from generic expectations of UV and IR divergences. Wether 
this also holds for the sub-leading in e parts of the calculation is a question which will be 
studied in this section. 

In general one should be looking to ways in which the powercounting argument pre- 
sented above could fail after integration. By the general argument just mentioned it is 
natural to suspect the UV and IR divergences of the integrals. The first arise from the 
momentum integration directly, while the second are related to the Feynman parameter 
integrals. To get a grip on the situation and also to be able to compare to results in the 
literature most of the concrete results in this section will be obtained at one loop. 

Another way of phrasing the problem under study is the question whether amplitudes 
are 'single-cut' constructible. As outlined above, the singularities of the integrand should 
be related to the branch cuts in the z plane. Any possible discrepancy arises from the 
integral over the arc at infinity. Since this is technically much easier, one would like to 
calculate cuts in four dimensions, not (4 — 2e) of them. This gives information on the 
leading poles of the e expansion, but disregards information about sub-leading terms. The 
question is if the information of the leading poles is enough to determine the full amplitude. 
In other words, the question is if boundary terms are missed if the limit e — > is taken 
before the z — > oo one. Based on known results [15] [7] it is easy to speculate that this 
possible problem is absent in supersymmetric theories with massless matter but is present 
in non-supersymmetric gauge theory. 

In the following we are in particular interested to find deviations from the general form 
of the tree level amplitude at large z as given in equation (2.23) as this would signal a key 
difference between integral and integrand. 

4.1 BCFW shift of the integrals at one loop 

In this subsection the powercounting argument will be applied to the (gluonic sector of the) 
one-loop quantum effective action as calculated in the background field method [36] [37]. 
The quantum fields will be put in background Feynman- 't Hooft gauge, while the tree level 
fields will be kept in AHK gauge. One direct motivation to study background field methods 
is that in this gauge in Yang-Mills theory it is easy to see supersymmetric cancellations 
at the one loop level directly. Moreover, it allows us to disregard the problems of the 
integration of unphysical poles in the momentum integration associated to the lightcone 
gauge, while keeping the analysis of tree level diagrams presented above intact. 
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The gauge invariant one-loop contributions to the effective action as calculated through 
the background field method read 



-^eff ^tree ^scalars ^fermions ^gluon j (4-1) 

with 

^scalar = log det;i Q (D^) (4.2) 

^chiral fermion = log det }_ 1 {D^D" + O^F^) (4.3) 

-Cgiuon = log det^ + E^F"") + log detj =0 (D^D") . (4.4) 



Here D„ is the covariant derivative for the external field A, a is the Lorentz generator in 
the spinor representation (a^ v = I [7^,7^]) and E is the Lorentz generator in the vector 
representation. The determinants are shorthand for the path integrals to be calculated 
which generate the Feynman graphs for the vertices. Explicit color-ordered Feynman rules 
can be found for instance in figure (8) in [38]. Explicit results for the three and four vertex 
as calculated through the background field method obtained with a different motivation 
can be found in [39] , [40] . 

In the Feynman- 't Hooft background field gauge the vertices of the quantum effective 
action are all proportional to integrals of the form 

k f d D l l^ 1 ... 



4 J (2*) D P(l + kl f ...(/ + EITi 1 ki) ' (45) 

where ki are the external momenta and there are k powers of I in the numerator. By 
convention the shifted momentum will be taken to enter the color ordered loop loop at the 
leg with momentum k\ and to exit at the leg with momentum &2 as the focus will primarily 
be on shifts of color- adjacent graphs. The generalization to more 'distance' between the 
shifted legs is trivial for the general argument to be presented below. D denotes the 
space-time dimension of the integrals, which will be taken to be D = 4 — 2e below 3 . In 
gauge theory k < n holds by simple power counting. Incidentally, note that the most UV 
divergent graph (k = n) consists of three point vertices only which is also the class of 
diagrams important to the large BCFW shift limit in Feynman- 't Hooft gauge at tree 
level. In the following some details will be needed about the structure of the integrations 
appearing at the one loop level. 

Momentum integrals 

The standard method of calculating the class of integrals in equation (4.5) proceeds through 
the introduction of Feynman parameters 

I* = Tin) / da n / y^4n — i , (4.6) 



3 We will opt to keep Lorentz indices in the numerator in D dimensions as this is more convenient in an 
off-shell formalism. 
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with 



J da n = m daij <5(^ Qi - 1) , (4.7) 

followed by a shift on the integration variable, 

p p + SP , (4.8) 

with 



n / i-1 \ 

«" = -E K (49) 

i=2 \j=l / 

which brings the integral into the form 

n = r( n) y da n y (2 ,/ ; /2+ ^ • 

In this expression 

A = 5>£> +( T) 2 , (4.11) 

i=l j=0 

where fco = is defined for convenience. At this point the momentum integral can be 
performed. This follows from the repeated use of 

_). Z 2 -^- , (4.12) 

which follows from spherical symmetry. This leads to the generic momentum integral 
formula 

dPi (-i)^nn-f-o <1 ..., 2i /iy-^ (413) 



(2tt) d (Z 2 + A) n (4vr) D / 2 2T(n) ' s V A 

where rfg 1 '"^ 2 ' is a sum over products of metrices with completely symmetrized indices. 
This leaves a generically complicated integration over the Feynman parameters. 

Feynman parameter integrals 

The integration over Feynman parameters is generically of the form 

I'lr ' (414) 

with / a polynomial of the Feynman parameters of maximal degree n and Sij is some sym- 
metric matrix which depends on external momentum invariants whose form is unimportant 
at the moment. There is a special sub-class of these integrals which can be integrated easily, 

/k k 
da n J] (a,) ft = r(fc + E-ft) II ^ + &)) , (4-15) 
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which holds as long as the arguments of the gamma functions all have positive real part, 



sft(l + A)>0 Vt. (4.16) 

At this point one can use the general approach of Mellin-Barnes integration (see e.g. [41] 
and references therein) to reduce the starting Feynman parameter integral to the special 
sub-class. The general formula needed to make this work is 

1 1 1 [ +io ° . X 



(X + Y) x T(\)2Tri 



/-\-tOG X' W 
dw yx^n™+^)n-™) , (4.i7) 
-ioo * 



where the integral is along a contour from —ioo to -Moo which splits the two series of 
poles of the Gamma function. This contour will be taken to be a straight line. By encir- 
cling either series of poles and summing the resulting residue integrals one can see that 
the above formula is closely related to Newton's binomial formula. Combining multiple 
applications of equation (4.17) and the integration of (4.15) one can see that one can trade 
the generic integral over Feynman parameters of equation (4.14) for multiple Mellin-Barnes 
type integrals. 

One difficulty of the outlined approach is that in general the constraints of equation 
(4.16) may be violated for small values of e. One solution to this is to first continue e to a 
value for which a good straight line contour can be found. Graphically, this puts the series 
of poles of all the involved T functions to the left (for T(x) terms) and right (for T(—x) 
terms) of the chosen contour. Then e is tuned back to a small value. This procedure will 
pick up residue integrals where poles cross the contour which can be calculated separately. 
This particular procedure is conveniently implemented in publicly available code [42]. 

Example: massless box 

As the Feynman parameter integrals may not be as familiar as the momentum integrals, 
perhaps an example is in order. For this, consider the Feynman parameter integral appear- 
ing in the calculation of the massless box, 



1 

(saiQ3 + ta2CXiY 



B = I da A - — • ( 4 - 18 ) 



Following the outlined program above one arrives at the following Mellin-Barnes type 
integral 

1 1 r+ioo w 

B = ¥^T(^2e)i_- dw \t) nw + 2 + e)T(-w)T(-w-l-efT( W + lf . (4.19) 

To make this integral well-defined a contour needs to be specified. It is not too hard too 
see a straight line contour with %t(w) = — \ for e = — 1 will neatly divide the poles of the 
6 Gamma functions. Continuing e to zero then picks up a residue integral at w = — 1 — e, 



1 \ s f 1 



/+ioo 
*. R ■ (4 ' 20) 
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for some function f of e and a remaining Mellin-Barnes integral over a remaining integrand 
denoted by R. The function / has a (4r) type pole typical of an IR divergence. Now 
suppose a BCFW shift of the box is considered of the type 

s ->■ s t^t + z', (4.21) 

In this case it is seen that the integral scales one order of z' better than the integrand. 

4.2 One loop momentum integrals contributing to large BCFW shift 

After this brief general discussion of the integrals appearing in the calculation of the one 
loop quantum effective action the next step is to discuss which one-loop diagrams contribute 
to the large BCFW limit of complete one loop amplitudes. For this one considers the one- 
loop vertices with hard lines entering on adjacent positions calculated using color ordered 
Feynman rules. In the following the focus will be on momentum integration, leaving the 
Feynman parameters for remarks at the end of the section. 

A simple class of diagrams is if the hard line only 'glances' the loop (that is, it connects 
to the loop through one of the four- vertices with two background fields). The one- loop 
vertex is of order (z°) as it is either proportional to a metric or is anti-symmetric. From 
the tree level one can simply use the tree level analysis to show that those diagrams where 
both shifted legs end on the loop vertex directly will not contribute new types of term to 
the generic large z-form of equation (2.23). A new class is formed by diagrams dressed by 
one order z° tree level hard leg as obtained from equation (2.19). For terms proportional 
to the metric this does not change the analysis, but for the anti-symmetric terms new order 
z° possibilities arise. It can be checked that this does not influence the 4 point amplitude, 
but it is an issue for higher points. 

This leaves diagrams where there is at least one loop edge along the hard line which 
requires more care. First consider the effects of the momentum integration. One can see 
from equations (4.9) and (4.11) that in general both the shift <5 M and the function A are 
linear in z, 

81* -> ^ - zazq" (4.22) 

n i—l 

A -> A - 2za 2 J2 a iC^2lj +S ) " 1 ■ ( 4 - 23 ) 

i=l j=o 

There is one loophole here where the two shifted external wave functions connect to the 
loop directly and form a triangle shaped diagram. The involved diagram topologies are 
illustrated in figure 3. In this case the two external momenta appearing in A are both 
orthogonal to q and hence A is in this case independent of z. Note that the same class of 
diagrams also contains those diagrams which are singular for the AHK gauge on the tree 
level parts as a sub-class. 

After momentum integration the ra-point vertices in the quantum effective action at one 
loop where the shifted legs connect to different vertices on the loop integral are proportional 
to 

[n/2] 
j=0 
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(a) (b) 

Figure 3: diagrams corresponding to a triangle shaped loophole to the scaling of equation (4.23) 

for some field dependent coefficients c& n . Here [n/2] is the greatest integer for which 
[n/2] < n. This allows us to classify the possible troublesome terms in D = 4 — 2e, which 
are at: 

n k z — order 
2 2 z 1 - 6 



2 z°~ e 

3 2 z 



o- e > ( 4 - 25 ) 



3 3 

4 4 z°- £ 

while the three particle vertices have a maximal divergence degree of z. Using equation (4.1) 
and taking into account the possible tree level contributions this table can be translated 
into diagrams again. The origin of the diagrams from the two different types of vertices 
in the determinants (either the D 2 or the F type) will be stressed as this allows an easy 
way of comparing theories with different matter content. As a complication, the tree level 
AHK gauge propagator has a potential z° piece, see equation (2.19). This is proportional 
to two g's though so this can only contribute in exceptional cases. In the figures below 
which correspond to the suspect diagrams this exceptional part of the AHK propagator 
along a hard line is indicated by a crossed out gluon propagator. 

The listed Feynman graphs are calculated in terms of integrals in appendix B, up to 
terms which vanish in the large z limit as analyzed through the momentum integrations. By 
tracing through the terms several results follow. The leading term in z is given by diagram 
(a) on the left hand side of figure (3) whose value is calculated in equation (B.2). This 
three particle diagram neatly parallels the diagram dominating at tree level. Structurally 
this term yields a contribution to the large z behavior of 

(4.26) 

Here fi^v is antisymmetric and the functions fi have expansion in e of the form 

fi ~ 1 + 0(e) , (4.27) 
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(/) (.«) (h) 

Figure 5: Diagrams connecting the shifted legs directly to the vertex 



from the Feynman parameter integrals. A - divergence in /2 arises from the Feynman 
parameter integral. 

The form of the exact large shift terms is found by contracting a tree level three vertex 
onto the one- loop three vertex. By the analysis of the gauge singularities there, the leading 
term in z arises by contracting two tree level currents with a metric and multiplying this 
with 

+ O(z ) , (4.28) 

where K4 is the total momentum of all particles connecting to the current J a {K4). The 
sub-leading terms arise from a variety of diagrams. It should be noted that the part of 
the above result proportional to k^^k^^ is highly reminiscent of the analysis of the BCFW 
shift of two gluons in the bosonic string in [4]. There it was shown that this particular 



A(z) ~ (J P (K 3 



kz^k 



3,1/ 



A:2 
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contribution is absent in the superstring. For the calculation here this can be checked as 
well: in the supersymmetric case this term is absent. 

4.3 One scalar loop 

For the purposes of this article attention will now be further restricted to the terms which 
involve a scalar loop only. This is done since these terms generate so-called rational terms 
at the one loop level. With this motivation BCFW shifts of these contributions have been 
studied as such in the literature as recalled in the introduction. This corresponds for the 
leading term written above to setting /q = 0, as can be verified straightforwardly from 
equation (B.2). This leaves the analysis of the sub-leading terms, written in equations 
(B.3) - (B.10), suitably restricted to the scalar contributions. 

In a class of its own is the contribution of the remaining triangle integral, diagram b, 
given in (B.3). This is of order z° from the outset. Furthermore, it is easy to see that the 
results follow the same pattern as the term just analyzed. As such it is again proportional 
to the metric from the UV divergent part, proportional to k^^k^ v for the UV and IR finite 
part or antisymmetric in \xv from the IR divergent part. The UV and IR finite part arises 
from the scalar loop and entails a metric contraction between the soft legs. 

This leaves those integrals for which there is a non-trivial dependence on a shifted mo- 
mentum in the loop integral. These diagrams can only contribute if they are UV divergent 
by the reasoning above. Note that this severely restricts the class of diagrams all the way 
down to four field terms. This is interesting in itself as it is known that the rational terms 
of the Yang-Mills amplitudes can be determined if the four field amplitude is known [8] [9] . 

Since one of the main interests in this section is to find deviations from the established 
pattern of (2.23) one can disregard any contributions proportional to a metric contraction 
of the shifted legs. In addition, terms antisymmetric in the shifted legs can be dropped as 
well. Another way of enforcing this is to study shifts of adjacent like helicity particles. As 

recalled in the appendix, the external wave functions for the (+- 1-) and ( ) shifts scales 

as - and z respectively and are to leading order the same. 

One quickly discovers all the terms of equations (B.5) - (B.10) are proportional to the 
same pre-factor for the like helicity shifts, 

44A$ P * ~(D - 2)e^T(e) ( ^"^ C) ) (z qK,f ^r,^) + O Q (4.29) 
= F pa = Fe^^p . (4.30) 

Here two overall factors F and F pcr have been defined for notational convenience and K4 
is the momentum on one of the two soft legs. This result holds regardless of any on-shell 
conditions on the external legs other than the shifted ones. What remains is to add the six 
diagrams to determine the correct pre-factor. Since one can quickly check in the four point 
case that the above equation would give a non-zero contribution to a shift of the four point 
all helicity equal amplitude it is natural to suspect that, when summed, all diagrammatic 
contributions add up to a function proportional to e to counter the UV pole. In fact, one 
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arrives at the following list of results, 



A c -»• -IF A d -»• i(3-2e)(l-e)F 

A e -> i(e - 1)F A/ -> |(3 - 2e)(2 - 2e)F 

A" -+ -(3-2e)F A h -> (2)F 

in 4 — 2e dimensions. When added, this yields 

confirming the earlier suspicion. A useful numerological trick to check consistency during 
the calculation is to track powers of 3 appearing as coefficients for the (-) poles: as 
these poles should cancel coefficients which do not multiply a three should cancel amongst 
themselves. This involves three instead of six diagrams. 

In summary, at leading (z 1 ) and sub- leading order (z°) at one (scalar) loop there are 
only very few diagrams contributing as analyzed from the momentum integration. One 
class is given by gluing the triangle loop to a tree level gluon current, while the other 
consists of two tree level gluon currents glued to the four vertex calculated above. All 
these terms are manifestly finite. The leading order term arises from one diagram only. 

Example: residues of helicity equal amplitudes 

The above can be used to compare to known results on BCFW shifts of single loop gluon 
amplitudes and in particular those in [7]. Here we will restrict to the helicity equal am- 
plitudes. For definiteness and without loss of generality we pick '+' helicity. The above 
results involve the tree level gluon current, which can be calculated in a variety of ways. 
For our purposes the all-plus gluon current as derived in [30] will suffice, 

J °°' 3+ ' 4+ --" + )~ W (34)...W ' (4 ' 32) 
In this expression q& is the dotted spinor associated to the light-like vector and K is 
the momentum in the off-shell gluon. Combining this with the leading order triangle loop 
derived above, one immediately obtains 

. . . ~ (fci - k 2 + zq) p J p + = , (4.33) 

for any number of external legs. Hence for an adjacent shift of like helicity particles on 
the like helicity amplitude the scaling is of order (z°) and originates from the dressed four 
vertices discussed above. Also diagram (b) on the right hand side of figure (3) whose value 
is calculated in equation (B.3) does not contribute as a metric contraction between two 
soft all-plus currents vanishes. This leaves diagrams which involve equation (4.31). 

For definiteness, consider a BCFW shift for legs 1 and 2 for which the shift vector is 

Qaa = la^a ■ (4.34) 

Therefore the gauge reference spinor in the current has to be set to 

C = 2° . (4.35) 
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The external wave functions of the shifted legs behave under the above shift as 

, &1 Q +0 (1\ (4. 36 ) 



(£1) (£(1 + 22) z(&) \z 2 

(£2) (£2) (£2) + U ^ ^ ' (4 - 3?j 

for some arbitrary gauge spinor £ which may not be taken to be equal to 2„. From the 
above the following form of the large z shifts follows, 

n— 1 



A(z) ~ ( g |^|^ J^(3, . . . , j)>(j + 1, . . . , n) j + O 0) . (4.38) 



It is easy to see that in this expression the £ dependence of the shifted wave functions drops 
out. In the case of four points there is only one term, 

«->~'OT +0 (i)- 



#)~-f|SV»f'] ■ (4-40) 



which is equivalent to 

4f^ F ( - 

[23] [14] 

This is indeed the known residue at infinity of the adjacent shift of the four point all plus 
amplitude at one loop, up to and including all the dependence on e. This follow because 
the amplitude is proportional to the following integral 

dD ( l )(v 2 ) 2 T?r, — r^Ti — r — r^ni — r^? » ( 4 - 41 ) 
P(l + fci) 2 (Z + k x + k 2 ) 2 (l - k 3 ) 2 

where n 2 is the invariant length of the part of the dimensionally regulated loop momentum 
which is orthogonal to the four dimensions of choice. This is indeed proportional to F in the 
limit, including all e dependence up to the subleading term in the factor of (D — 2) = 2 — 2e. 
The latter discrepancy can be explained from e scalar terms usually ignored in amplitude 
calculations at one loop. In the four point case the residue at infinity coincides with the 
complete amplitude. 

For arbitrary multiplicity at leading order in e, the expression in equation (4.38) reads 

A ^~* 1<2 3 }<34}..W (£g WtjTU) J + °UJ +0(e) ■ 

(4.42) 

Now for a fixed value of k, say a, and a fixed value of I, say b, there are only some specific 
terms which contribute to this expression 

(...k = a,i = b . ..) = (, w I , ^gM^oy + i)^)^], 

1 ' V(23)(34)...(n2)M^ (2j)(j + 1,2) 
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This sum can be performed using 

^(2i)(i + l,2) (2a) (62) " 1 " j 

Now the total contribution to the residue is the sum over all possible a and b such that 
a < b. This reads 



This can indeed be checked to be the residue at infinity as derived from the known expres- 
sion at leading order in e. 

4.4 Remarks 

It is quite plausible that the above derived large z boundary terms are complete in the sense 
that for pure Yang-Mills these are the only boundary terms to be found. At leading order 
in z this indeed follows from the background field analysis in this section. The sub-leading 
order requires still more effort. However, all terms from the diagrams in figures 3, 4 
and 5. which potentially contribute at order (z°) and which are not accounted for by the 
scalar contribution would contribute - IR poles to one of the purely rational amplitudes in 
pure Yang-Mills at large BCFW shift. Since it is known these are absent in the full answer 
they must cancel from other sources. The divergences arising from Feynman parameter 
integrals have been left unevaluated until this point. As will be briefly indicated below, 
these divergences require care. 

For supersymmetric theories it is clear that the leading pole in z at one loop at least 
has the same form as at tree level (being proportional to a metric). In fact, from the way 
this pole arises from the special three point vertex diagram it is expected that this can be 
proven fairly easily to all orders in perturbation theory. This follows since the form three 
particle coupling is very much restricted in Yang-Mills especially when supersymmetry is 
taken into account. Again, the sub- leading orders may be very complicated and the real 
problem is to show that there are no sources of z dependence from all the other diagrams 
which interfere with this. For maximally supersymmetric theories knowledge of one good 
shift is enough to show that good super-shifts exist, see [34]. 

Although a particular gauge setup was employed in deriving them there should be 
a more gauge invariant formulation of the same. After all, the full amplitudes are gauge 
invariant. For this one can note that additional BCFW shifts on other legs may be employed 
to derive auxiliary recursion relations for the order z contributions. For the rational terms 
this is the strategy advocated in [8] [9]. What the above adds is a deeper understanding 
why these recursion relations work once the four point amplitude is known. 

Note that the analysis gives a neat interpretation where the difference between cuts in 
four and higher dimensions is for pure Yang-Mills: if e — ?• is taken first for cuts in four 
dimensions the rational parts are found at the residue at z = oo. If the residue at z = oo is 
calculated first, the rational amplitudes are at the branch cuts. This quantifies the order 
of limits problem mentioned above. 
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Improved shifts arising from IR divergences 

From the massless box example one can see how terms which have been counted as sub- 
leading in z up until now may be elevated. The point is that if possible divergences are 
encountered when deriving the Mellin-Barnes representation, there will be an enhancement 
of the z scaling. See for instance equation (4.20) and the discussion below which illustrates 
this point neatly. On the other hand, for the massless box type integral this only arises if 
the exponent i in 



is bigger than 2 and there are no redeeming powers of Qj around to dampen the divergence. 
Moreover from inspecting the poles passing the contour this effect can only elevate a naively 
order z~ % suppressed integral to -. Hence this mechanism only contributes at higher order 
in z if there is another, explicit power of z floating around in the integral which can be 
used to generate a contribution and the integral is sufficiently IR infinite. This cannot 
happen for the integrals appearing in the scalar loop. For the complete background field 
calculation for the gluon in the loop note that for the leading term in z this can at most 
contribute a metric as the two possible sources of z in the numerator of the integrand 
are the momenta at legs 1 and 2. To yield a z = (j) x z 2 contribution these have to be 
contracted somewhere else in the diagram and not to each other. This leaves the metric 
contraction between the shifted legs. Intriguingly, this term is proportional to \ whose 
numerator on general grounds is expected to be related to the tree amplitude again and it 
is encouraging that the leading poles indeed conform to this expectation. 

A full discussion of IR divergences will be deferred to future work. For inspiration see 
[43] and references therein. 

5. Towards BCFW on-shell recursion relations at loop level 

Although this article mainly deals with the large BCFW shift behavior of the integrand 
and its integrals one of the prime objectives of this line of research is to obtain similarly 
useful on-shell recursion formulae at loop level as exist at tree level. In this section some 
(very) preliminary developments in this direction are presented. The main question is 
what the singularities of the integrand and the integrals correspond to in terms of lower 
loop or lower point amplitudes. Since branch cut singularities in the integral generically 
correspond to a cut propagator in the integrand these questions are essentially the same, 
up to the discrepancy quantified in the previous section. This will be ignored here as it is 
mostly relevant when discussing the dimension in which the cuts are calculated. 

From the good large z behavior for color- adjacent shifts one can write a prototype 
recursion relation for any amplitude. 




(4.46) 





where the superscript I in A® denotes the loop order. Here particles one and two have 
been shifted, which are assumed to be color adjacent. The first term arises from poles in 



tree parts of the amplitudes and are just a natural generalization of tree level BCFW, i.e. 
Al) _ sr^ ST'ST' _(^> • • • ^ Pr)AlJ_ k+2 (Pr, k +!,...,!) 

^tree poles" 2^ 2^ 2^ (h + h\ 2 ' [ ' 

r,h(r) i=0 fc=3 (K 2 + . . . KfcJ 

where the first sum is over all different mass levels r and over all polarization states at that 
level, denoted h{r). The new term arises from single cuts of the amplitude and reads 



n) , (5.3) 



c,r,h(r) 

where c sums over the possible single cuts and I is on-shell, 

f -'-«sTi- (5 - 4 » 

The question therefore is if there is an interpretation of the quantity B^^ 1 ^ in terms of 
lower loop, higher point amplitudes. A natural guess would be 

Bfc 2 1) (iJ.-U,3... J n) 4+2 } (i,r,-U,3...,n) . (5.5) 

Note however that from the diagrams it is easy to see that B is finite while the term on 
the right hand side would be evaluated at a collinear singularity. More cunning is therefore 
required to resolve this issue. The simplest suggestion would be to simply subtract of the 
collinear pole from the lower loop amplitude in some regularized fashion as it is basically 
the product of two amplitudes again, but this would leave a soft pole to deal with. See [21] 
for some more suggestions how to define this 'forward limit' properly in a very much related 
context. Any of these suggestions immediately leads to a concrete recursion relation. Note 
that the authors of [22] sidestep the issue by working only with those MHV diagrams which 
do not have tadpoles. This is another way of defining the quantity B to any loop order 
(see also [44]) and would express the tree amplitude in terms of tree amplitudes again. 



Figure 6: Diagram topologies contributing to the four point scattering amplitude at one loop in 
<f> theory 

Some further inspiration can be taken from eft 4 theory. There the four particle scattering 
amplitude for instance has two diagrammatic contributions at one loop as depicted in figure 
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6. Sums over different particle orderings will be suppressed. BCFW shifting particles on 
two adjacent legs gives a (z°) behavior of the integrand, exactly as at tree level. This can 
be sidestepped by introducing more particles which modify the UV behavior of the theory 
[35] but what is important here is the pole of the integrand which does get a shift. Cutting 
the associated line gives a tree structure, obviously. However, this is not a 6 particle tree 
level amplitude: what is missing for this is the Feynman graph which has the cut legs 
attaching to the same vertex on the tree. This tree diagram can be added and subtracted 
schematically as, 

(Ai)({ Pi }) = J ( g g 6D{P ± At ({ft}, I, -l) - F ({ft}, I, -l) . (5.6) 

When employed within the recursion approach the integral over the second term would 
vanish in dimensional regularization as it is a tadpole type integral. For Yang-Mills this 
suggests to sum only those diagrams which do not have the cut legs attached to the same 
vertex, which naively coincides with the prescription in [22]. 

Nested shifts 

One obvious suggestion to sidestep the forward limit issue is to use an auxiliary BCFW 
shift on two adjacent legs on the residue of the single cut. If the residue is gauge invariant 
up to terms which integrate to zero, one can use the same argument as above to argue 
that the single cut residue again scales as a shift of a tree level Yang-Mills amplitude to 
all loop orders. Now at one loop for instance this reduces the single cut residue to either a 
more familiar double-cut residue or contributions which involve B type terms with less legs 
dressed by tree amplitudes. At higher loops part of the terms obtained will be double-cuts 
while others will be two single cuts of separate loops. Iterating shifts should eventually yield 
an expression of the amplitude in terms of integrals over tree level amplitudes evaluated 
at shifted momenta. 

Schematically this reads at one loop, 

B <°> 2 ( U -1 2,3 ...,»)- E E ^•-"~^ < f r ^" g --) , ,5.7) 

r,h(r) fc=3 \ > 

where the tildes denote the second BCFW shift, taken on particles n and n — 1. The 
quantity C is either an A or a B depending on the location of the I legs: if both of them 
are on a certain C it is a -B-type, if one of them is both C's are an A and if none of them 
are on a certain C it is an ^4-type. Now equation (5.1) can be used in reverse to express 
the -B-type terms on the right hand side again in terms of lower (maximally n — 1) point 
one loop amplitudes. This expresses a one loop scattering amplitude in terms of scattering 
amplitudes only with either a lower loop or a lower leg order. 

The resulting expressions are rather messy and at higher loops the involved structures 
for nested cuts get more and more intricate. It will be very interesting to see if more 
efficient formulas than this can be made to work, especially when combined with maximal 
supersymmetry. Note that in deriving the formulas above several issues with order of limits 
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problems between the two shifts and a potential problem with assumed gauge invariance 
have been neglected: these most certainly deserve further study. 

6. Discussion and conclusions 

In this article a first step was made towards on-shell recursion at loop level. It was proven 
the Yang-Mills integrand can be reconstructed from its singularities in full analogy with 
tree amplitudes by studying large BCFW shifts of the integrands. Clearly, the shifts are 
probing a form of universal behavior in Yang-Mills theory. The main obstacle to apply 
on-shell recursion relations is to find a proper interpretation of the singularities of the 
integrand in terms of integrands for lower loop or lower point amplitudes. The preliminary 
analysis above for double shifts leaves much to be desired and clearly much more interesting 
work can and needs to be done here. 

It should also be very interesting to extend the analysis of the difference between 
integrand and integral made explicit above at one loop. Especially higher loops would be 
most welcome. That this should be possible follows from the link to UV and IR singularities 
of the integrand. For UV singularities at one loop a full discussion was given above. The 
next step here is to understand the shift behavior of the Feynman parameter integrals at 
one loop in more technical detail, which should enable the natural extension of the analysis 
to two loops and beyond. It is clear that the number of Feynman graphs to be calculated is 
going to be limited, especially compared to usual two-loop computations. This should yield 
already some interesting new results for the general form of two-loop amplitudes. Note that 
the analysis of large BCFW shifts gives information on complete amplitudes which is in a 
real sense complimentary to the usual unitarity cuts. It would also be interesting to see if 
this information can already be used to explicitly calculate amplitudes. For instance, the 
large z behavior obtained above could be used in principle to constrain the rational terms 
to higher (all) orders in e. 

From the results above it is easy to speculate that the Einstein gravity amplitude 
integrand will obey the same large- z scaling behavior as at the tree level (i.e. the 'square' 
of that of Yang-Mills theory). This is certainly plausible both from the string theory point 
of view as well as tree and loop level experience. This leads one to suspect for instance that 
loop level relations for integrands like those conjectured in [45] between gauge theory and 
gravity can be proven just as the tree level version of these relations [46] can be proven from 
non-adjacent BCFW shifts [47] (see also [48]). In this respect it is encouraging that non- 
adjacent shifts of particles on the same color trace for Yang-Mills amplitudes are suppressed 
by Q) for all cases considered in this article by simple absence of the leading pole diagrams. 
This implies the existence of loop level 'bonus' relations. Moreover, the results on shifts 
obtained in this article is independent of the color structure of the unshifted parts of the 
diagrams: they hold in particular for non-planar amplitudes. 

It will be interesting to study the application of the ideas of on-shell recursion at loop 
level to maximally supersymmetric Yang-Mills theory in four dimensions as it opens a 
window to prove conjectured properties of the integrals and integrands such as dual con- 
formal invariance and its one-loop breaking (see [49] and references therein and thereto). 
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Also, there should be a close link to on-shell recursion in twistor spaces and the conjec- 
tured Grassmanian structure of the leading singularities of the amplitudes in maximally 
supersymmetric Yang-Mills. It is exceptionally natural to transform the above recursion 
formulae and its natural generalizations to super twistor and super ambi-twistor space in 
the light of the results in [5] and [6]. This is in a real sense the same intuition which 
motivated the present paper (see subsection 3.1): on-shell recursion should be related to a 
form of the CFT bootstrap equations also at loop level. Twistor spaces are after all closely 
connected to the conformal nature of M = 4 super- Yang-Mills. Actually this is not too far 
off the idea that also at weak coupling maximally supersymmetric Yang-Mills should be 
related to a string theory [1]. This can perhaps be made more precise starting from the 
tree level results in [50]. 

Finally, it should be noted that there no inherent obstructions to the loop order of 
the amplitudes under study for the techniques developed here. This opens the window to 
NNLO calculations in QCD for instance an admittedly very small fraction. Certainly first 
the one-loop structure needs to be under full control and cross-checked to known results 
before any serious attack can be made on experimentally relevant calculations beyond the 
one loop level. However, as methods beyond one loop for standard model calculations 
are currently limited to delicately applied brute force even a slightly opened window is a 
tempting prospect. 
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A. BCFW shifts of fermions in D dimensions 

In this appendix the analysis of shifts of gluons for the Yang-Mills integrand in the main 
text is extended to shifts of external fermions in D dimensions. In addition this appendix 
spells out some conventions about quantum numbers and the BCFW shift. The starting 
point is 



for two particles labelled 1 and 2. As noted in the main text, the following vectors span a 
four dimensional subspace, 




k 2 - qz 



(A.l) 
(A.2) 



Span ( 



q,q,h,k 2 ) 



(A.3) 



- 32 - 



of a D-dimensional space. In other words, there is a choice of normalization of q and a 
choice of Lorentz frame such that 

g=-^(0 ) l,*,0 } ...) (A.4) 

g=-L(0,l,-i,0,...) (A.5) 

fc 1 = A(i )0 ,...,l) (A.6) 

fc 2 = A(i >0 ,...,-1) . (A.7) 

Moreover, a choice of ordering of the entries has been used to make the formulae more 
readable. These equations are for instance the starting point of the higher dimensional 
analysis of BCFW shifts in [26]. From the explicit expressions it is easy to see that q and 
q are helicity eigenstates of momenta k\ and k2, where helicity is defined w.r.t. the four 
dimensional subspace, i.e. 

Ri = 1 \ M T , (A.8) 

where e^ 4 ) is the totally antisymmetric symbol in the subspace spanned by the vectors in 
(A. 3) and E is the generator of rotations in D dimensions. Note that e^ 4 ) depends non- 
trivially on the momenta in general. Also, this is the helicity operator for k2 while it is 
minus the helicity operator for k\. 

It is very natural to take the momentum of the other leg as the choice of gauge for each 
of the legs. The polarizations of the gluon are then given by q and q and can be split up 
into either nontrivial (±) or trivial (transverse) under the above helicity generator. Under 
the BCFW shift these obey 

9 1 = at = <? -> 9 

9 2 = Q -> Q ~ k k i 

T T 

9 -> 9 

This concludes the analysis of the external wave functions for gluons. What is needed 
below in addition to this are the polarization spinors. These are solutions to the Dirac 
equation for either k\ or k2 with definite eigenvalues under the generator (A.8). In the 
special four dimensional frame these can be constructed using the four dimensional spinor 
helicity method. However, in higher dimensions in general it will pay to be slightly more 
general, this will be done below. 

The main observation needed to study the spinor polarizations is that one can construct 
the set 

{ q ,q,k 1 ,k2,gTgf} , (A. 10) 

with all vectors light-like to form a basis of M. D with D even. Here the transverse polariza- 
tion vectors have been grouped into (normalized) conjugate pairs, 

gf.gj = S ij i,j E {!,..., D- 4} . (All) 
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rf 


= 73^ 
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= 




1 T 

= V2 9i 7 ^ 


r r+i 


1 -T 

= V29i 7^ 



In passing we note the relation of the set just constructed to the choice of complex structure 
on MP which can be seen most easily from equation (A. 4). Using the basis of (A. 10) one 
can study representations of the 7 matrix algebra, 

{7^,74 = 2^,, (A.12) 

by projecting. Define for this the operators 



(A.13) 



so that the Clifford algebra reduces to D/2 copies of the (normalized) fermionic harmonic 
oscillator, 

{r+ rj} = l . (A.14) 

The standard representation theory of this algebra can be written neatly in terms of the 
half-integer eigenvalues under the 'rotation' generator 

Ri = r-r+ - \ . (A.15) 

Let the vector of eigenvalues be h. All spinors are annihilated by half of the generators, 
while the other half of the generators flips one eigenvalue, 



rf\h) = 6 hi ^_\h') , (A.16) 



where h' has the i-th eigenvalue (labelled hi) inverted. The chirality of the spinor is the 
product of the signs of the eigenvalue labels. Of course, the first eigenvalue label on the 
spinors corresponds to the spinor being a solution to the Dirac equation for k\ for /io = +5 
or hi for h§ = — \- The constructed spinors form a basis of the space of all spinors. 

With this construction in hand, it is now straightforward to study the BCFW shift 
considered in this article. The shift concerns only the first two quantum numbers of the 
spinors as expected and leads to the following shifted spinors 

II I \ -> II I ) 

I 2 ' 2 ' ' ' 7 I 2 ' 2 ' ' ' 7 

I I 1 \ v. |1 1 \ _ z I _ 1 1 \ 

I 2' 2' -- 7 ~~ M 2' 2' ■ ■ 7 k I 2'2>"7 /A i7\ 

i 1 1 v . I _ 1 _1 v , z |1 1 \ • \ A - L <) 

I 2' 2i ■ ■ 7 — f I 2' 2' - • 7 "r fe I 2' 2' ' - 7 

III \ I _ 1 1 \ 
I 2 ' 2 ' - - 7 ^ I 2 ' 2 ' • • 7 

It is not too hard to verify that expressed in terms of four dimensional spinor helicity 
this is the usual BCFW shift. The four dimensions of course correspond to the special four 
dimensional space spanned by k\, k2, q, q and the first two quantum number labels on the 
spinors. The four dimensional chirality /helicity of the spinors for instance is the product of 
the signs of these two labels. The remaining labels can be interpreted as the R-symmetry 
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labels in this four dimensional decomposition. The four dimensional interpretation can 
hence be captured in the following notation 

il>j(ki) -> i>7(h) 

^{k 2 ) -> ^j{k 2 ) 

where the R-symmetry indices should be interpreted as a (D — 4)/2 dimensional vector 
with half-integer entries. This is the notation used below. 

For completeness the shifts of the conjugate spinors can be derived from the complete- 
ness relation, 

t„(*i + *<iY = ■ (A.i9) 

In particular, 



h=± 



^f(h) -> iPi(h) 

^7(fci)-^7(*i)-f V>7(fc2) 



V>7(fc)^i (to) + 



(A.20) 



^/ (to) -> V>/ (to) 

is obtained through this. The conjugate spinors are the solution to the conjugate Dirac 
equations. The above analysis is an application of techniques explored in [25]. Indeed, the 
above can easily be turned into a completely covariant treatment by replacing 

k -> y/ki ■ k 2 (A.21) 

everywhere and phrasing the basis in terms of its inner products only. 

BCFW shifts involving one or two fermionic legs 

Now the stage is set to extend the reasoning in the main text to shifts of one gluon and one 
fermion and two fermion legs in generic minimally coupled (renormalizable) gauge theories. 
For shifts involving a fermion a choice of gauge is not necessary as the naive powercounting 
can be used to derive the large z behavior, as shown below. For fermions the analysis below 
makes the argument touched upon in [52] explicit. 

Two shifted gluons 

As shown in the main text for two gluons the large z behavior of the amplitude follows 
from 



A = el (% x ) e v 2 (k 2 ) = 4 (z) e v 2 (z) Urj^z + B^ u + (J) ) . (A.22) 



Here e± and e 2 are the polarization vectors of the shifted adjacent legs, the hats on the 
momenta denote shifted momenta and f\ is an arbitrary function of the unshifted momenta 
and polarization vectors. The matrix B is antisymmetric. This Ansatz follows from the 
above argument about the class of dominating diagrams. Combining the large z-behavior 
of the polarization vectors with the Ansatz yields table 1. 
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One shifted gluon, one shifted fermion 

In general diagrams containing external fermion legs have a definite ordering of the fermion 
lines which has to be picked. This determines the external fermion wave functions to be 
given by either ifi or its conjugate ip. The quantum numbers of the fermions can be selected 
too, see the appendix for a construction. 

In this case the large z behavior of the amplitude follows from 



A = et{ k x M^u (k 2 )= el (z) r£ + K a <f + q a K* + hw° 



7m (1 + K'wf] +o(^\\u (z) , (A.23) 



where u(k 2 ) is the (in general z-dependent) polarization spinor of the fermion, ip is an 
arbitrary (conjugate) spinor, f\ is an arbitrary function and K, K and K are arbitrary 
vectors. 

In general the hard line connecting the external particles contain scalar fields, multiple 
fermion fields and gluonic contributions. The leading contribution is again a combination 
of very simple three point vertices. Combinations of hard gluon and scalar lines coupling to 
this glue have been analyzed above. Note that possible vertices with the Yukawa couplings 
between scalars and fermions are sub-leading in the large z limit as they are momentum 
independent. There is at least one fermion attaching to the hard line since this connects to 
the external particle. This fermion vertex lowers the overall z count by one. The vertex to 
which the shifted gluon leg attaches directly on the hard line contains one gamma matrix 
in the complete gamma matrix trace, whose structure can be denoted as 

tl®lu{i)l^...{i))u{k 2 ) . (A.24) 

In this expression the fa indices on the gamma matrices attach to other parts of the diagram 
which are left undetermined as they will not influence the final result. Since (^) 2 = this 
reduces to 

'^(l + ^))u(fe)f]]J , (A.25) 

by repeated use of the gamma matrix algebra. Additional fermions along the gluonic hard 
line will lead to an additional suppression of - for diagrammatic reasons. Assembling the 
hard gluon line contribution with the above coupling to the fermionic part of the hard line 
now leads to the equation (A.23). The large z-behavior of the amplitudes which follows 
from the ansatz will be listed below in table 2. 

Two shifted fermions 

For the case of two shifted fermions the large z behavior of the amplitude follows from the 
following 2 particle current, 



A = u (k^jMu (k 2 ) =u(z) (k 1 ^ + - z K 2 + LLX\ u (z) , (A.26) 
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-1 


+ 1 


-1 


-1 
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Table 2: leading asymptotic power in z K of the adjacent BCFW shift of two particles in a tree 
amplitude in D dimensions for all possible polarizations 

which follows directly from power-counting the involved diagrams. The z-independent but 
otherwise arbitrary functions K\ and K 2 can only contain an even number of gamma 
matrices. The first arises from a fermionic hard line while the second contains one gluon 
exchange diagram. 

The large z behavior which follows from the above analysis of shifted fermions is 
captured in table 2. The entries containing one fermion and one gluon in this table were 
discussed in [52]. In the table the fermion polarizations are listed as defined in a special 
set of 4 dimensions and split into helicity and R-symmetry labels as explained more fully 
above. The difference between the entries for g T and g T2 on one axis and fermions on the 
other is that 

sjyy^O while sJVV^O, (A.27) 

while the two different entries for the fermion shifts correspond to either the same R- 
symmetry labels (1) or different ones (0). Note that the fermion entries in this table were 
derived by straightforward powercounting, while only the sub-leading behavior of the gluon 
shifts required more care in the form of the above AHK gauge powercounting. 

As a consistency check at tree level one can use the supersymmetric Ward identities to 
check that if the last two columns of table 2 are known, the other entries may be derived. 
The power of the supersymmetric Ward identity is that it is independent of the coupling 
constants in the theory, so these relations must hold for amplitudes at any loop level in 
Yang-Mills theory. 

B. One loop diagrams contributing to the BCFW shift at leading order 

In this appendix the explicit results for the one loop diagram contributions as obtained 
through the color ordered background field method used in the setup in the main text are 
collected. These have been calculated with the color ordered Feynman rules as can be found 
in figure (8) in [38]. Legs one and two will be taken to be shifted. Sub-leading terms in 
the BCFW shift parameter z will be discarded as indicated. Subscripts on square brackets 
indicate how many powers of the background field type vertex coupling two quantum fields 
to the background field field strength have been used. Schematically this reads 



{a,F^[A]a u y o [...], 



(B.l) 



The diagrams themselves are listed in figures 3, 4 and 5. 



Triangle vertices with two shifted legs 

This class of diagrams is illustrated in figure 3. The trivalent triangle vertex evaluates to 

d D l 
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(2vr)^ l2{l- kl f{l + k 2 Y 



{(D-2) [iM 2l + k 2-h) P \ 



+4 [(l/j.rjpv + IvVpp ~ lpV/j,u)(fa ■ k 2 )j 
The triangle vertex with two shifted legs and two external legs evaluates to 
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Here the definition 



was used. 



k = fa + k 2 



(B.3) 
(B.4) 



Bubble and triangle vertex contributions with one or more tree level propaga- 
tors 

This class of diagrams is illustrated in figure 4. It consists of bubbles and triangles con- 
nected to the loop vertices through the order z° part of the lightcone gauge propagator. 
Many terms vanish because this part of the propagator is proportional to q fJ, q p , see equation 
(2.19). This implies only one or two three vertices with a shifted leg can be sewn onto the 
loop vertices. 

The contributions of a trivalent triangle vertex with one added lightcone gauge part 
evaluates to 
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The bubble vertex with one 'lightcone leg' contributes 
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with a | symmetry factor included here and in every following bubble. The bubble vertex 
with two 'lightcone legs' contributes 



e l e 2 A p,vpu — e l e 2 



J P L^^^^_(\_ l m -^\ q a {2l + k 4 ) a q^{2l + k i )p 



Four field bubble, triangle and box one-loop vertices 

This class of diagrams is illustrated in figure 5. The four field bubble evaluates to 
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The two (mirrored) four field triangles evaluate to 

fPfZAte) =fVf u [ - 

1 2 ^ vpu 1 2 J (2n) D P(l + k 1 ) 2 (l-k 4 ) 2 
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Finally, the four field box evaluates to 

/dPl 1 

[27:rP(i + k l ni+k 1 +k 2 ni-k A f • • ^4 

+4 [(^/p^ - lj p rj^)(h ■ k 4 )} 2 + (D-2) [l p lj p l a } ) + O (^) (B.10) 

at leading order in z. To simplify the second term it was noted that under the loop integral 

l^lvVpcT = IplaVpu + O Qj (B.ll) 

holds. Furthermore, 

(ki ■ h) = -(k 2 -k A ) + (z°) (B.12) 

was used. 
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